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Abstract 

We measure the influence of different time-scales on the dynamics of financial 
market data. This is obtained by decomposing financial time series into 
simple oscillations associated with distinct time-scales. We propose two new 
time-varying measures: 1) an amplitude scaling exponent and 2) an entropy¬ 
like measure. We apply these measures to intraday, 30-second sampled prices 
of various stock indices. Our results reveal intraday trends where different 
time-horizons contribute with variable relative amplitudes over the course of 
the trading day. Our findings indicate that the time series we analysed have a 
non-stationary multifractal nature with predominantly persistent behaviour 
at the middle of the trading session and anti-persistent behaviour at the open 
and close. We demonstrate that these deviations are statistically significant 
and robust. 

Keywords: Scaling exponent, Entropy, Hilbert-Huang transform. Market 
efficiency. 


1. Introduction 

Financial markets are complex and dynamic systems that generate non- 
stationary, non-linear and noisy time series. These time series do not consist 
of a single driving force, but of several components that are superimposed 
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onto each other in a hierarchical form. The oscillating components are gen¬ 
erated by many participants with diverse interests and investment horizons 
that interact and create effects that seem to repeat themselves cyclically with 
different characteristic periods [1]. These effects can be retrieved and quan- 
tihed by looking at the statistical properties of hnancial data at different 
time-scales. 

A classical approach to model hnancial markets has been the efficient 
market hypothesis (EMH) [2], which states that hnancial markets are efficient 
if they rehect all available information, thus arbitrage conditions are quickly 
eliminated. According to this theory, stock prices are unpredictable. The 
weak form of the EMH admits a rapid price adjustment process [2]. However, 
in practice, prices tend not to adjust to new information so rapidly, taking a 
certain amount of time. During this time investors can take actions to exploit 
temporary prohtable opportunities arising from new information [3, 4]. 

According to the fractal market hypothesis (EMH) [5], hnancial partici¬ 
pants are heterogeneous and market stability exists if there are investors with 
different time-horizons that create liquidity. Market participants have invest¬ 
ment horizons that vary from seconds to years (market makers, noise traders, 
hedge funds). They treat the arriving information differently and affect the 
price dynamics in various ways depending on their trading time-scales [6]. 
The existence of investors with different investment horizons induces a stable 
market, markets become unstable when one horizon becomes dominant since 
liquidity dries up. 

The EMH predicts that critical events are connected to dominant in¬ 
vestment horizons and is usually associated to self-similarity, fractality or 
multifractality. Self-similarity in hnancial data was hrst studied by Mandel¬ 
brot [7, 8], and has been found to be present across hnancial markets with 
complex properties that are signihcantly related to economic and hnancial 
characteristics of the markets [9, 10, 11, 12, 13]. Self-similarity is related 
to the occurrence of similar patterns at different time-scales. In this sense, 
probabilistic properties of self-similar processes remain invariant when the 
process is viewed at different time-scales [7, 14]. A stochastic process X{t) 
is statistically self-similar, with scaling exponent 0 < 77 < 1, if for any real 
a > 0 it follows the scaling law: 

X{at) = a^X{t) t G M, (1) 

where the equality (=) is in probability distribution [15]. 
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An example of self-similar process is fractional Brownian motion (fBm), 
a Gaussian process with stationary increments characterized by a positive 
scaling exponent 0 < H < 1 [16]. When 0 < H < 1/2, the increments 
of fBm show negative autocorrelation. The case 1/2 < H < 1 corresponds 
to fBm with increment process exhibiting long range dependence, i.e., the 
autocorrelation of the increment process is positive. When i/ = |, the fBm 
reduces to Brownian motion (BM), a process with independent increments 
[15]. 

It has been observed that each moment of the distribution of hnancial 
returns varies as a power law of the time horizon with a different H exponent 
[14]. This phenomenon is called multiscaling and reflects the occurrence 
of different dynamics at different time-scales, it can be attributed to the 
heterogeneity of market participants. Time-dependent scaling behaviour has 
also been observed in hnancial time series [17, 18], the local variations of 
roughness can be described by allowing the H exponent to vary with time 
[19]. 

Entropy measures have also been used to measure the complexity of h- 
nancial times series [20, 21, 22]. A low value of entropy indicates the presence 
of more predictable patterns that are therefore associated with periods of h- 
nancial inefficiency. Conversely, when the time series exhibit more irregular 
and less predictable patterns, the uncertainty level is higher and such periods 
are described by larger values of entropy. 

In this paper, we propose two new time-dependent measures based on the 
scaling properties of hnancial time series. These novel measures are obtained 
via the Hilbert-Huang transform [23], a methodology that as a hrst step, 
decomposes the analysed time series into several oscillatory modes by means 
of empirical mode decomposition (EMD). Secondly, the Hilbert transform 
is applied to these oscillations to obtain time varying attributes. The time- 
dependent scaling properties of hnancial time series are associated with the 
relative weights of the amplitudes at characteristic frequencies. 

The hrst measure that we introduce in this paper is a scaling exponent 
that quantihes the relative hierarchical variations of the amplitudes of the 
components with respect to their associated time-scales; the second is an 
entropic measure quantifying the dispersion of the amplitudes of the compo¬ 
nents. 

The remaining of this paper is organized as follows. In Section 2, we 
introduce the methodology used to identify the oscillating components of 
the data, the Hilbert-Huang transform. The proposed scaling measure with 
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some application to self-similar process is described in Section 3. In Section 
4, we propose a comparative entropy-like measure. In Section 5, we apply 
the proposed measures to intraday financial data. Finally, in Section 6, we 
provide conclusions and future perspectives. 

2. Hilbert-Huang Transform 

The Hilbert-Huang transform (HHT) [23] is a technique used to analyse 
non-linear and non-stationary time series. It was originally designed to study 
water wave evolution, but it has proven to be a useful tool for other complex 
signals, including hnancial time series [24, 25, 26, 27, 28]. 

The HHT consists of two steps: namely, empirical mode decomposition 
(EMD) and Hilbert transform (HT). The EMD decomposes the time series 
into a set of narrow-band intrinsic mode functions (IMFs) and the Hilbert 
transformation of these IMFs provides local frequency and amplitude at¬ 
tributes. 

The EMD is a fully adaptive decomposition that in contrast to Fourier 
and wavelet transforms, does not require any a priori basis system [29]. Fur¬ 
thermore, it can be used to analyse non-stationary time series, it assumes 
that any time series consists of super-imposed oscillations. The purpose of 
the method is to identify these oscillations with scales dehned by the lo¬ 
cal maxima and minima of the data itself. Hence, given a time series x{t), 
t = 1,2, ...,T, the EMD process decomposes it into a hnite number of IMFs 
denoted as Cfc(t), k = l,...,n and a residue function, r(t). The number of 
IMFs, n, is approximately of the order of log 2 (T) [23]. The IMFs are com¬ 
ponents oscillating around zero which are obtained through a sifting process 
that makes use of local extrema to separate oscillations starting with the 
highest frequency. At the end of the sifting process, the time series x{t) can 
be expressed as: 


x{t) = '^Ck{t)+r{t). (2) 

k=i 

The residue function, r(t), is the non-oscillating drift of the data. For more 
details about this decomposition refer to [23]. The EMD was proposed as 
a way to pre-process time series before applying the Hilbert transform. It 
generates components of the time series whose Hilbert transform can lead to 
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physically meaningful definitions of instantaneous amplitude and frequency. 
The Hilbert transformation of each function Ck is dehned as: 


OO 

cl{t) = - [ dr, (3) 

TT 7 t - T 

— OO 

where the integral has a singular point at r = t and it is dehned as a Cauchy 
principal value. Each function Ck and its Hilbert transformation, Cfc, produce 
a complex function Zk, dehned by Zk{t) = Ck{t) + Ck{t) with amplitude 

afc(t) = ^/ck{ty + Ckity, (4) 


and phase 


Ok{t) = tan 


1 Ck{t) 
Cfc(t)’ 


( 5 ) 


The instantaneous 
respect to time 


frequency is dehned as the derivative of the phase with 


UJk{t) 


d0k{t) 

dt 


( 6 ) 


The instantaneous amplitude results in a smooth function, an ‘envelope’, 
that takes the overall shape of the time series taking the maxima and minima 
but never crossing the time series itself. This amplitude function represents 
the combined oscillations of all frequencies involved in the component but 
removes all the frequency information. 


3. Time-dependent Scaling Exponent 

The proposed time-dependent scaling exponent, denoted as is con¬ 

structed by observing the way local amplitudes ak{t) (Equation 4) change 
with respect to the local periods Tk(t) = (Equation 6) for all k = 

1,2,...,u. 

We hrst applied the method to fBm and observed that the amplitude 
function obtained through the HHT follows a power-law behaviour with re¬ 
spect to the instantaneous period: 

ak{tf (X (7) 

where the exponent H*[t) describes the local scaling properties of the IMF 
amplitudes, and is comparable in magnitude to the self-similar exponent H. 
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Logr^(3421) LogTj5405) 


Figure 1: Illustration that for fBm the local amplitudes ak{t) and the periods Tfc(t) follow 
Equation 7: ak{t)^ oc Plots report instantaneous amplitude as a function 

of period for the following four randomly chosen times: t = 1326,2252,3421,5405. The 
simulated process is a fBm with self-similar exponent H = 0.6 and length T=10,000 points. 
The straight lines are best-fit linear regressions. 


In Figure 1, we report a particular instance of Equation 7 showing the 
linear ht between logafc(t) and logTfc(t) for four randomly chosen times of a 
fBm with self-similar exponent H = 0.6 and length T=10,000. The values 
of H*[t) reported in the plots are obtained from the slope of the regres¬ 
sion line. We observe that they are all consistently close to the self-similar 
value H = 0.6. For the chosen values of t, we calculated the goodness of 
the linear ht by estimating the coefficient of determination [30] (values 
of this coefficient range from 0 to 1, with 1 indicating a perfect ht between 
the data and the linear model). Results for the following four randomly 
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chosen times: t = 1326,2252,3421,5405 are as follows: i?^(1326) = 0.99, 
i?2(2252) = 0.90, i?2(3421) = 0.98, i?2(5405) = 0.99, indicating therefore 
that for those instants of time the data are well represented by the log-linear 
model of Equation 7. Similar linear scaling results are obtained across all 
times, but the scaling exponent is different at each time step, making H*{t) 
a time-dependent estimator. A value of H*[t) > 0.5 is obtained when around 
time t, the amplitude of long cycles is larger than in a pure random walk. 
This can be interpreted as a persistent behaviour in the amplitudes of the 
process, meaning that in a neighbourhood of time t the process is in a cycle 
indistinguishable from a trend. On the contrary, values of H*{t) < 0.5 rep¬ 
resent a rougher and more chaotic behaviour around time t. These processes 
are composed of oscillations with more similar amplitudes across time-scales 
than in Brownian motion, creating a complex and uncertain behaviour. In 
this case, high frequency components are more active and their contribution 
to the total variance is more signihcant than in a random walk process. 

3.1. Extended Simulation of Self-similar and Long Memory Proeesses 

To test the power-law relation of Equation 7, we extended the simula¬ 
tion set of fBm and we considered other two different self-similar processes, 
namely a-stable Levy motion (SLM) [15] and autoregressive fractionally in¬ 
tegrated moving average (ARFIMA) processes [31]. 

For each stochastic process, we simulated m = 1,000 paths of length 
T = 10,000 points^. We estimated H*{t) and calculated the time-dependent 
sample mean over the number of simulations, i.e., we calculated {H*{t)) = 

hTZiHtit). 

We also estimated the sample mean of H*[f) over time and over the 
number of simulations, {{H*)) = L {H*{t)). The standard deviation of 
is calculated as: 

= v'eL. E™, ww - irnyf /(T-m- 1). 

For each process, the average over time and over the number of simu¬ 
lations of the coefficient of determination is also estimated and denoted as 
((A^)). We compared the estimated H*{t) with the generalized Hurst expo¬ 
nent [33] with q = 1 and denoted as Hq. 


^The length of the LFSM is set to T = 2^® — 6000, following the algorithm proposed 
in [32] 
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Figure 2: a) Illustration that for £Bm the scaling exponent H*(t) is on average close to 
the self-similar exponent H. The plot reports the sample mean of denoted as 

{H*{t)) and computed over m = 1,000 simulations of fBm with self-similar exponent 
H = 0.1,0.2,..., 0.9 (bottom to top) and length T = 10, 000 points, 
b) Sample mean of H*[t) over time and over number of simulations, denoted as {{H*)), 
standard deviation of H*{t) and sample mean of the coefficient of determination, 

The values of GHE(l) denote the generalized Hurst exponent with q = 1. 


• Fractional Brownian motion. Stochastic processes with scaling ex¬ 
ponent varying from H = 0.1,0.2, ...,0.9 were simulated. All the sim¬ 
ulations were done using the Matlab® wavelet toolbox. Results for 
different values of H are reported in Figure 2(a). We observe that 
{H*{t)) consistently varies around the input self-similar value of H. 
Some disturbances can be observed at the beginning and end of this 
time series due to the boundary effects of the EMD [34]. 

In Table 2(b), we report {{H*)) and the standard deviation of 
We observe a good agreement with the self-similar parameter H, but 
large values for sidn* that could be attributed to the local characteris¬ 
tics of H*{t). We notice that for fBm with scaling exponent H < 0.3, 
the ((R^)) coefficients are small, indicating significant deviations from 
the scaling law of Equation 7. We obtained consistent results when 
comparing the estimated H*{t) with the generalized Hurst exponent. 

Moreover, we also considered fBm paths with shorter length, T = 1, 000 
and T = 500 points, we do not report these results, but we noted that 
the longer the time series, the better the estimation of the scaling ex¬ 
ponent H. Likewise, the standard deviation and the goodness of the 




















fit improve with the length of the time series. However, all results are 
consistent with the ones reported here for length T = 10,000. 


• a-stable Levy motion (SLM). This process is a generalization of 
the Brownian motion to the a-stable distribution with 0 < a < 2. The 
case 0 = 2 corresponds to Brownian motion. The SLM is a l/a-self- 
similar process {H = - ) with stationary and independent increments 
[15]. Similarly, the extension of the fBm to the a-stable distribution 
is the linear fractional stable motion (LFSM), a self-similar process 
that exhibits both heavy tails and serial correlation. The increments 
of LFSM have long range dependence when H > l/a and negative 
dependence when H < l/a [15]. Though the a-stable Levy motion 
is self-similar with H G [1/2, oo), the self-similar parameter of LFSM 
varies in (0,1). For more details about these processes refer to [15]. 

We generated SLM processes using the toolbox provided by [32], sample 
paths are of length T = 10, 384, with parameters for the generation 
m = 128 and M = 6000, making m{M -|- T) to be a power of 2, see 
[32] for more details. We considered the case H = 1/a for values of 
H = 0.5,0.55, ...,0.95. 

The time-dependent sample mean over the number of simulations is 
displayed in Figure 3(a). We observe a noisier estimator than the one 
obtained for fBm. In Table 3(b), we report {{H*[t))), noticing a fair 
approximation to the self-similar parameter H, with better results for 
processes with H <0.7. The means of the coefficient of determination 
suggest that the scaling relation of Equation 7 is indeed satished. We 
compared the proposed estimator with the generalized Hurst exponent 
with q = I, obtaining consistent results, i.e., Hq = 1/a [35]. 

• ARFIMA. We tested the log-linear relationship of Equation 7 in 
ARFIMA(p, d, q) processes with Gaussian innovations, p,q eN, autore¬ 
gressive and moving average coefficients respectively [31]. This model 
is an extension to ARIMA(p, d, (;) models, allowing the differencing 
exponent d to take fractional values, —1/2 < d < 1/2. The correspon¬ 
dence between the two parameters H and d is given by id = d -|- 1/2. 
The interval 0 < d < 1/2 of long range dependence corresponds to 
1/2 < id < 1 [36]. We considered the simple case of ARFIMA(0,d,0) 
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Figure 3: a) Illustration that for SLM the scaling exponent H*{t) is on average close 
to the value H = -. The plots report the sample mean of H*{t), denoted as {H*{t)) 
and computed over m = 1,000 simulations of SLM with self-similar exponent H = 
0.5,0.55,..., 0.95 (bottom to top) and length T = 10,000 points. 

b) Sample mean of H*[t) over time and over number of simulations, denoted as {{H*)), 
standard deviation of H*{t) and sample mean of the coefficient of determination, 


with fractional order d = —0.4, —0.3,..., 0.4 and length T = 10,000. 
We calculated H*[t) for the integrated ARFIMA time series. 

From Figure 4(a), we observe that the estimator {H*) is a good ap¬ 
proximation of the exponent H. In Table 4(b), we report {{H*)), the 
standard deviation of H*[t) and the sample mean of the coefficient of 
determination, {{B?)). Similar to the fBm case, the estimation is more 
accurate for FI > 0.3 and the coefhcient of determination is closer to 1. 

From the analysis of these three different processes, we observe that our 
proposed method produces a fair approximation to the self-similar parameter 
H. However, let us remark that this scaling exponent is not intended as an 
alternative method to estimate H, which can instead be obtained with more 
reliable tools [11, 37, 38]. The aim of this method is instead to compute the 
time-dependent amplitude contribution of the prevalent fluctuations present 
in a time series, distinguishing between periods when high or low frequencies 
are contributing more or less than what could be expected for Brownian 
motion. 
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Figure 4: a) Illustration that for ARFIMA(0,1,0) the scaling exponent H*{t) is on average 
close to the value H = d + 0.5. The plots report the sample mean of denoted as 

{H*{t)) and computed over m = 1,000 simulations of ARFIMA(0,d,0) with self-similar 
exponent H = 0.1, 0.2,..., 0.9 (bottom to top) and length T = 10, 000 points, 
b) Sample mean of H*[t) over time and over number of simulations, denoted as {{H*)), 
standard deviation of H*{t) and sample mean of the coefficient of determination, 


4. Time-dependent Complexity Measure 

We define a time-dependent Shannon entropy-like measure based on the 
square of the amplitude of the IMFs. This measure provides a time-varying 
quantihcation of complexity that offers an alternative to the scaling exponent 
to measure the strength of cycles present in hnancial time series. Making use 
of the functions a^, described in Equation 4, we dehne a time-scale relative 
distribution of amplitudes as: 

Ph*) = (8) 

E “») 

k=l 

where n is the number of IMFs excluding the residue. Similarly to Shannon 
entropy [39], we dehne the time-dependent complexity measure as: 

n 

C*{t) = - Inpfc(t). (9) 

k=l 

Equation 9 provides a measure of the distribution of amplitudes between 
the oscillating components. If the total amplitude at time t is concentrated 
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in one oscillation mode, we observe a low complexity value, implying that 
around time t the process is following a prevalent trend. On the contrary, if 
at time t all the oscillation modes have similar amplitudes, we obtain a large 
complexity value that indicates a more erratic and unpredictable behaviour. 

Thus, C*{t) provides a time-varying estimation of disorder and adapts 
closely to our visual perception of complexity. Moreover, Equation 9 offers a 
more general measure of uncertainty than the variance since the latter mea¬ 
sures the dispersion around the mean, while C*{t) measures the dispersion 
of energy around the different IMFs. Similar to an entropy measure, the 
value of the proposed complexity at time t varies between zero, if one IMFs 
dominates the energy of the process, and log(n) if the energy is uniformly 
distributed between the n IMFs. 

The choice of weights equal to the square of the amplitudes in Equation 8 
is arbitrary, although it is in agreement with other measures of entropy that 
have been dehned, for example in [40]. We tested alternative choices, such 
as the linear weight ak{t), obtaining analogous results to the ones reported 
here. 

Let us note that, although not independent, the two measures convey 
different information. The estimator C*{t) is an information quantiher of 
uncertainty, it is obtained from the distribution of the amplitudes regardless 
of their time-scales and it only quantihes the homogeneity of the components. 
On the other hand, the scaling exponent, measures the change in 

the amplitudes across time-scales, testing the scaling law in Equation 7. In 
this respect, H*{t) is a more restrictive measure that assumes a log-linear 
relationship between amplitudes and periods. 

5. Time-dependent Scaling in Financial Markets 

We applied the proposed measures to intraday prices of four stock indices: 
(1) S&P 500 (USA), (2) IPC (Mexico), (3) Nikkei 225 (Japan) and (4) XU 100 
(Turkey). We intentionally chose two hnancial markets that are classihed as 
developed (USA and Japan) and two emerging markets (Mexico and Turkey) 
with the additional feature that the Japanese and Turkish stock exchanges 
have two trading sessions separated by a lunch break. 

The data set, obtained from Bloomberg, consists of prices recorded at 
30-second intervals. It covers a period of 5 months, from January 15*^, 2014 
to June 16*^, 2014. The number of days and the number of data points 
for every trading day depend on the opening hours of each stock exchange. 
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(a) s&p 500 (b) ipc 



(c) Nikkei 225 (d) XU 100 


Figure 5: 30-second sampled log-prices for different stock indices for the period January 
15*'*, 2014 to June 16*'*, 2014. (a) S&P 500, (b) IPC, (c) Nikkei 225 and (d) XU 100. 


Logarithmic prices of each hnancial stock index are plotted in Fignre 5. Table 
1 shows the number of days and the length of each time series. 


Country 

Index 

No. of Days 

Length 

USA 

S&P 500 

105 

81,900 

Japan 

Nikkei 225 

104 

62,400 

Mexico 

IPC 

101 

78,780 

Turkey 

XU 100 

106 

78,440 


Table 1: Number of days and length of each financial time series. 

The time evolution of H*(t) over the 5-month period for the four hnan¬ 
cial indices is shown in Figure 6 (dark-blue line). We note that H*[t) has 
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large intraday variations that make difficult to identify any trend over longer 
periods. For this reason, we report a moving average version of H*{t) de¬ 
noted as H*{t) (light-blue line in the same Figure). More specihcally, H*{t) 
is calculated from the relation, afc(t)^ oc where akit) and fkit) are 

the averages over a rolling window of the size of a trading day. The dashed 
red line in this Figure indicates the value H = 0.5. 



0.2'- ^^^^ - 0.2'-^^^- "I- 

15/01/14 13/02/14 17/03/14 15/04/14 15/05/14 16/06/14 15/01/14 14/02/14 14/03/14 14/04/14 15/05/14 16/06/14 


(a) s&p 500 (b) ipc 



15/01/14 13/02/14 14/03/14 15/04/14 16/05/14 16/06/14 15/01/14 13/02/14 14/03/14 14/04/14 15/05/14 16/06/14 

(c) Nikkei 225 (d) XU 100 

Figure 6: Time-dependent scaling exponent for different stock indices for the period Jan¬ 
uary 15*^, 2014 to June 16*^, 2014. The scaling exponent H*{t) is depicted by a dark-blue 
line. The light-blue line represents a rolling average over the length of a trading 

day. The red line indicates the value H = 0.5. (a) S&P 500, (b) IPC, (c) Nikkei 225 and 
(d) XU 100. 

By comparing the values of H*{t) and H*{t) for the four different stock 
indices, we observe that the S&P 500 index is the one closest to the value 
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H = 0.5, that assuming a Gaussian distribution will imply BM behaviour, 
see Figure 6(a). In this Figure, H*(t) fluctuates around H = 0.5 exposing 
only some brief departures from it. For instance, we can detect a period 
around February 2014 where the scaling parameter results in signihcantly 
larger values. In this period of time, the S&P 500 index was indeed in a rising 
rally, see Figure 5(a). Therefore, this suggests that the identihed persistent 
behaviour could be attributed to the recognition of a long time-scale cycle 
with larger amplitudes than in the case of pure random walk. 

In Figure 6(c), we report the scaling dynamics for the Nikkei 225 index. 
We observe that H*{t) has values constantly above 0.5, specially at the end 
of the analysed period. It should be noted that this market has lunch breaks 
that affect the intraday values of H*{t). 

For the IPC index the values of H*{t) and H*{t) are consistently closer 
to H = 0.6 than to the BM value H = 0.5, see Figure 6(b). This suggests 
that IPC returns show intervals where the amplitude displays a persistent 
behaviour. Similarly, the Turkish scaling exponents take values larger than 
H = 0.5, see Figure 6(d). 

We tested the validity of Equation 7 when applied to hnancial data by 
computing for every time t the coefficient of determination R^[t). The mean 
over the whole period is reported in the second column of Table 2. We also 
considered the three cases: H*{t) < 0.45, a window around 0.45 < H*{t) < 
0.55 and H*[t) > 0.55. We observe that the goodness-of-ht is generally 
better for H*[t) > 0.5, the interesting case when hnancial data show trending 
behaviour, see Table 2. 


Index 

{R^)au 

{R?) H*<0A5 

(-^^)o.45<if<0.55 

{R'^)h*>0.55 

S&P 500 

0.8753 

0.825 

0.8716 

0.8916 

IPC 

0.8812 

0.7971 

0.8703 

0.8915 

NIKKEI 225 

0.8072 

0.7345 

0.7829 

0.8198 

XUlOO 

0.9196 

0.7987 

0.8737 

0.9209 


Table 2: Average goodness-of-fit coefficient {R^) for the amplitude versus period log-linear 
model for different financial indices. First, the average is calculated for all the times t, 
then it is calculated separately for those times where H*{t) < 0.45, 0.45 < H*{t) < 0.55 
and H*{t) > 0.55. 

For a comparative analysis, we calculated the complexity measure C*{t) 
described by Equation 9. The obtained values for each hnancial stock index 
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(a) s&p 500 (b) ipc 



(c) Nikkei 225 (d) XU 100 


Figure 7: Time-dependent complexity measure, C*(t), for different stock market indices 
for the period January 15*^, 2014 to June 16*^, 2014. (a) S&P 500, (b) IPC, (c) Nikkei 
225 and (d) XU 100. 


are illustrated in Figure 7. The complexity values for S&P 500 index, Figure 
7(a), are overall the largest among the four indices. 

The IPC index shows an increasing evolution of snggesting a more 

uniform distribution of amplitudes at the beginning of 2014, Fignre 7(b). On 
the contrary, the Nikkei 225 index presents a decreasing measure of complex¬ 
ity, indicating a period of higher complexity at the beginning of the sample 
period, Fignre 7(c). Finally, the XUlOO index presents alternate intervals 
of high and low complexity, displaying regnlarly large values in the last two 
months of the sample period. This higher randomness is also visible from the 
scaling exponent that displays relative lower values of Figure 7(d). 
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Overall, the functions H*[t) and C*{t) vary in opposite directions. For 
each market, the correlation between these two measures is negative with 
values pskP = —0.21, pipc = —0.23, pNikkei = —0.28, pxu = —0.15. We note 
that these correlations values are small, indicating a weak linear dependence 
between these variables. This is to be expected as the underlying measures 
are associated with rather different properties. 

5.1. Intraday Analysis of Scaling Patterns 

We investigated the intraday patterns by separating the paths of H*(t) 
and C'*(t) into daily windows. Taking for example the time series H*{t) 
for the S&P 500 index. Figure 6(a), we separated this time series into the 
n = 105 days that compose the data set, see Table 1. In Figure 8(a), we 
display these daily time series (one day on top of the other). The colour 
bar represents the value of H*{t). This graphical representation allows us to 
compare trading sessions and identify patterns at specihc times of the day. 

We estimated the statistical mean of H*{t) across the days, resulting in 
an average value for each time t of the trading session. This average, denoted 
as {H*(yt ))describes the regular behaviour of H*{t) on a trading session, 
see Figure 8(b). 

In order to validate that the observed dynamics of {H*it))^^y^ are sta¬ 
tistically signihcant, we compared these dynamics with scaling exponents, 
obtained from several realizations of Brownian motions of length 
equal to the analysed hnancial time series, see Table 1. The time series of 
were fragmented into n windows of equal length. The mean over 
these n days is denoted as {H*bm{'^)) days- The pink band reported in Figure 
8(b) corresponds to the 5*^ and 95*^ percentile of the empirical distribution 
of {H*bm{'^)) days computcd from 100 simulations. 

We compared the values of H*{t) obtained for each hnancial index with 
the {H*bm{'^)) days band. At each time t during the trading session, we esti¬ 
mated the relative fraction of H*[f) values that falls outside the pink band. 
In Figure 8(c), we report these results as a ratio of number of days outside the 
band divided by the total number of days. This ratio is labelled as likelihood. 
The colour bar of this hgure represents the value of the average scaling ex¬ 
ponent, i.e, {H*{t))days (value plotted in Figure 8(b)). From this Figure, we 
observe that across the day there are periods of time with very high empirical 
probability of observing values of the scaling exponents signihcantly different 
from the corresponding values extracted from pure Brownian motion. 
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The mean of the complexity measure, at each time t of the trading 

session is shown in Figure 8(d). Equally as with the H*{t) exponent, the 
mean of C*{t) is computed across all n days and it is denoted as 
The same daily analysis for the remaining three financial indices is reported 
in Figures 9, 10 and 11, respectively. Overall, from the intraday scaling and 
complexity measures we observed the following patterns: 

• For each of the selected financial indices, the daily average 
displays an inverted U-shaped form that reflects a more chaotic be¬ 
haviour at the beginning and at the end of the trading session. The 
opposite behaviour is observed for the complexity measure, 

which reveals a U-shaped form. 

• The S&P 500 index displays the largest values of H*(t) (a stronger am¬ 
plitude persistent behaviour) during the middle of the trading session. 
From Figure 8(a), we observe that most of the days present large values 
of H*(t) around midday. At this time, the average exponent {H*{t))^y^ 
reaches a value of 0.6, see Figure 8(b). These values are signihcantly 
different from what would be expected for Brownian motion with more 
than 80% of the observations outside the 5*^ and 95*^ percentile. Figure 
8(c). Consistently, the complexity measure reaches its minimum at the 
same time of the trading session, see Figure 8(d). 

• The Mexican stock exchange is characterized by some large scaling 
values at the middle of the day. However the most noticeable pattern 
is the large values just before the end of the trading session, see Figure 
9(a). The mean of the windowed values reaches a maximum of 0.65, 
creating an upswing shape at this time, see Figure 9(b). This could 
be associated with an increase of trades in the last few minutes of the 
trading session that creates a more drastic change in the amplitudes. 
This pattern is only present in the Mexican stock exchange. 

From Figure 9(c), we observe that some minutes before the closing of 
the market, more than 90% of the local scaling indices fall outside the 
5*^ and 95*^ percentile band for Brownian motion. The complexity 
measure also reflects a steep increase of disorder at the end of the 
trading session, see Figure 9(d). 

• The Japanese and Turkish stock exchanges display two regions of large 
values for the scaling exponent. These regions are separated by the 
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Figure 8: Intraday analysis for the S&P 500 Index. 

(a) Intraday dynamics of the scaling exponent, H*{t), as a function of day and time. The 
colour bar indicates the value of the scaling exponent 

(b) Mean of H* (t) over the 105 days, denoted as {H* The pink band corresponds to 

the 5*^ and 95*^ percentile of the distribution of computed from 100 simulations. 

(c) Likelihood of H* (t) to fall outside the 5*^ and 95*^ percentile band for Brownian motion 

(pink band of Figure (b)). The colour bar indicates the value shown in 

Figure (b). 

(d) Mean of the windowed complexity measure, denoted as {C*{t)) 


lunch break, see Figures 10(a) and 11(a) respectively. The mean of 
the scaling exponent reflects a quasi-double inverted U-shaped form 
that is associated with the opening and closing of the morning and 
afternoon sessions, Figures 10(b) and 11(b). It is worth noting that 
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Figure 9: Intraday analysis for IPC. Caption for sub-figures (a), (b), (c) and (d) same as 
Figure 8. 


the two trading sessions do not display exactly the same prohle. For 
the Japanese hnancial market, the inverted U-shaped form of the hrst 
trading session is slightly skewed to the right, in comparison with the 
more symmetric shape of the second trading session, see Figure 10(b). 

• For the Turkish stock exchange, we observe that the hrst part of the 
trading session presents larger values of More than 90% of the 

analysed days present local scaling exponents that surpass the value of 
0.6, see Figure 11(c). The dominance of one IMF amplitude in the hrst 
trading session is also rehected in the lower values of the complexity 
measure, which reaches the lowest value when compared to the other 
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Figure 10: Intraday analysis for Nikkei 225 index. Caption for sub-figures (a), (b), (c) 
and (d) same as Figure 8. The white vertical band in each sub-figure corresponds to the 
lunch break in this stock exchange. 


stock markets, see Figure 11(d). 

Overall the intraday patterns of H*[t) and C*{t) confirm the well known 
fact that activity on hnancial markets is not constant throughout the day. 
The uncovered patterns corroborate the hectic bny and sell activity affecting 
different markets at the opening and closing of trading sessions [41, 42], 
Smaller values of H*{t) (large values of C*{t)) imply a non-persistent and 
rougher behaviour that is reflected in higher volatility. The exposed daily 
patterns of and are in agreement with the results that 

document the existence of a distinct U-shaped pattern in market activity and 
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Figure 11: Intraday analysis for XU 100 index. Caption for sub-figures (a), (b), (c) and 
(d) same as Figure 8. The white vertical band in each sub-hgure corresponds to the lunch 
break in this stock exchange. 


volatility over a trading day, i.e., volatility is higher at the opening and closing 
of trading and low in the middle of the day, see for example [43, 44, 45]. 

By comparing the day-by-day complexity values across the four markets, 
we observe that S&P 500 index displays the largest values across all the 
trading session. Nikkei 225 is the second most complex, followed by IPC 
and lastly the XUlOO, which has the smaller complexity values. This is in 
agreement with the results reported for developed and emerging markets 
[11, 46, 22]. 
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6. Conclusion 

We studied the relative weight of the oscillating components present in 
intraday hnancial time series coupled with their characteristic time-scale. 
These components are extracted via the Hilbert-Huang transform. We have 
shown that the combination of EMD and its associated Hilbert spectral anal¬ 
ysis offers a powerful tool to uncover the time-dependent scaling patterns of 
intraday hnancial data. 

We proposed two new time-dependent measures: 1) an amplitude scaling 
exponent and 2) an entropy-like measure. Using these measures, we have 
been able to identify trends and intermittent behaviour in hnancial time se¬ 
ries. Our measures are non-parametric and they do not assume any a priori 
stochastic process. The scaling exponent only assumes the existence of a 
power-law relation between the instantaneous amplitudes and the instanta¬ 
neous periods that was empirically shown to be present. 

We applied the scaling and the entropy-like measures to the study of four 
hnancial markets, two developed (US and Japan) and two emerging markets 
(Mexico and Turkey). We contrasted and compared the decomposition of 
their hnancial indices. With the use of intraday data, we recognized some 
patterns and identihed periods of low and high complexity. 

Compared to the other studied indices, the S&P 500 index results the 
most complex market. The intraday analysis reveals a distinctive anti- 
persistent behaviour at the opening and closing of the trading session, con¬ 
trasting with the persistent behaviour at the middle of the session. Similar 
intraday results are obtained for the other stock indices. The variations ob¬ 
served in the scaling and complexity measures are well outside the 5*^ and 
95*^ percentile expected for Brownian motion, suggesting strong deviations 
from this model that could be attributed to the presence of long range de¬ 
pendence or/and heavy tails. 

With the proposed measures, we are able to describe the dynamics of 
hnancial time series whose regularity changes over time. Our results sug¬ 
gest that hnancial time series have dynamic scaling properties that could be 
attributed to the autocorrelation of the process, the presence of tails and 
the non-stationarity of the time series. Given the presence of this dynamic 
roughness, a time-varying scaling exponent could better rehect the scaling 
behaviour of hnancial data. 
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